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The flat, homogeneous, and isotropic universe with a massless scalar field is a paradigmatic model 
in Loop Quantum Cosmology. In spite of the prominent role that the model has played in the devel- 
opment of this branch of physics, there still remain some aspects of its quantization which deserve a 
more detailed discussion. These aspects include the kinematical resolution of the cosmological singu- 
larity, the precise relation between the solutions of the densitized and non-densitized versions of the 
quantum Hamiltonian constraint, the possibility of identifying superselection sectors which are as 
simple as possible, and a clear comprehension of the Wheeler-DeWitt (WDW) limit associated with 
the theory in those sectors. We propose an alternative operator to represent the Hamiltonian con- 
straint which is specially suitable to deal with all these issues in a detailed and satisfactory way. In 
particular, with our constraint operator, the singularity decouples in the kinematical Hilbert space 
and can be removed already at this level. Thanks to this fact, we can densitize the quantum Hamil- 
tonian constraint in a well-controlled manner. Besides, together with the physical observables, this 
constraint superselects simple sectors for the universe volume, with a discrete support contained in 
a single semiaxis of the real line and for which the basic functions that encode the information about 
the geometry possess optimal physical properties. Namely, they provide a no-boundary description 
around the cosmological singularity and admit a well-defined WDW limit in terms of standing waves. 
Both properties explain the presence of a generic quantum bounce replacing the classical singularity 
at a fundamental level, in contrast with previous studies where the bounce was proved in concrete 
regimes -focusing on states with a marked semiclassical behavior- or for a simplified model. 
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I. INTRODUCTION 

Loop Quantum Gravity (LQG) [H-Q is one of the most promising candidates for a theory of quantum 
gravity. Nonetheless, although this quantization program of General Relativity has been extensively 
developed, it has not been completed yet. In order to test the program and obtain physical predictions in 
simple situations of interest. Loop Quantum Cosmology (LQC) implements the quantization procedures 
of LQG in symmetry reduced models This field has suffered a major breakthrough in recent years 
thanks to the large number of studies carried out in models such as flat FRW cosmologies [5l-[Tl| , FRW 
universes with other topologies [l^ - [T^ or with a non- vanishing cosmological constant [H a^Sotropic 
Bianchi I models [TgI - IiqI . or even some inhomogeneous scenarios po| . 

In this work, we will focus our attention on the flat FRW model with a massless scalar field, which is 
the most prominent system studied so far in LQC. In fact, even though this is the simplest cosmological 
model with non-trivial dynamics, its polymeric quantization has already provided relevant results [5|-|7[, 
such as the validity of the classical dynamics for semiclassical states in LQC in the region far away 
from the classical singularity and, more importantly, the replacement of this big bang singularity by a 
quantum bounce. In addition, the robustness of these results has been carefully checked [s^^. Besides, 
those studies, together with other complementary analyses have served to rigorously establish 

the mathematical foundations of LQC. 

Nevertheless, despite the extensive analysis of this model performed up to date, there are still some 
issues which are not completely clear and need a more careful discussion in order to have a thorough 
understanding of the theory. In particular, the points that we have in mind involve, on the one hand, a 
rigorous densitization of the Hamiltonian constraint (with respect to the volume of the universe) at the 
quantum level -a step which is carried out in order to obtain a densitized Hamiltonian constraint that 
is easier to solve than the non- densitized one arising naturally from LQG-, and, on the other hand, the 
superselection of the kinematical Hilbert space in different sectors which are intended to be as simple as 
possible while possessing optimal physical properties. Among such properties, we are interested in the 
existence of a regime with a well-defined WDW limit and where general features of the states provide a 
basis to explain the occurrence of a quantum bounce. 

With the aim of investigating these questions while keeping a rigorous mathematical control of all the 
steps in the discussion, we will adopt a new prescription when promoting the Hamiltonian constraint 
of the cosmological model to a symmetric operator. Actually, our quantum Hamiltonian constraint 
results more suitable than previous proposals inasmuch as it clarifies the issues stated above. One of 
the noticeable features of our prescription is the resolution of the cosmological singularity already at the 
kinematical level, in the sense that our Hamiltonian constraint allows us to decouple the zero-volume 
state from the rest of states in the kinematical Hilbert space, so that the kernel of the volume operator 
can be removed from the quantum theory (in a certain sense, this implements ideas from Bojowald 
[2]|). Once the singularity has been eliminated, we are able to formally establish a bijection between the 
solutions of the densitized and non-densitized versions of the quantum Hamiltonian constraint, an aspect 
which was not properly considered in previous analyses. Another important feature of our Hamiltonian 
constraint is that it decouples states with opposite orientations of the triad. Owing to this fact, we 
can identify in our theory superselection sectors which are simpler than those obtained in Refs. 0, 
since they are supported in discrete semilattices contained in a single semiaxis of the real line, instead 
of the whole line. Therefore, in our theory we can restrict the study to (e.g.) the positive semiaxis in 
a natural way, without the need to introduce any particular boundary condition or appeal to the parity 
symmetry discussed in Ref. (this symmetry is straightforward to impose in our case). The simplicity 
of our superselection sectors allows us to gain physical intuition and improve our knowledge about the 
physical consequences of the loop quantization. On the one hand, the basic elements which encode the 
information about the geometry, namely the generalized eigenfunctions of the gravitational part of the 
constraint (defined as a self-adjoint operator), are each explicitly determined by a single piece of initial 
data. They "arise" in a single section of minimum non-zero volume, without crossing the singularity to 
the sector of opposite triad orientation and without satisfying any kind of boundary condition. In this 
sense, the self-adjointness of the gravitational part of the constraint leads to a no-boundary description. 
On the other hand, whereas this part of the constraint is an operator with non-degenerate spectrum, its 
analog in the WDW theory is two-fold degenerate. This property, together with the fact that it is a real 
difference operator in the volume representation, implies that its eigenfunctions have a WDW limit with 
the form of an exact standing wave, i.e., each eigenfunction converges to a specific linear combination of 
the two analog WDW elements (one of them contracting in volume and the other one expanding) where 
they both contribute with equal amplitudes. Since, owing to the no-boundary property, the outgoing flux 
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associated with the expanding component must become incoming flux corresponding to the contracting 
component, and vice versa, a quantum bounce must happen. In this way, we will be able to explain 
the existence of a quantum bounce, showing that it is a direct consequence of the quantum geometry 
underlying LQC and hence a fundamental feature of the theory. 

The paper is organized as follows. In Sec. |ll] we review the classical framework and the quantum 
kinematics. We choose a new, suitable symmetric ordering for the Hamiltonian constraint and densitize 
it in Sec. IIIII The gravitational part of the constraint is analyzed in detail in Sec. IIVI We then impose 
the constraint in Sec. |Vl obtaining the physical Hilbert space. In Sec. IVII we compare different aspects 
of our proposal and of previous ones, showing the advantages of our approach. Finally, in Sec. IVIII we 
discuss the main conclusions of our analysis. 



II. THE FLAT, HOMOGENEOUS, AND ISOTROPIC MODEL 



We briefly summarize here the classical formulation of the model and its kinematical treatment in 
LQC. For more details, see Refs. [1, 0]- 



A. The Classical Approach 



In the flat model under study, the spatial hypersurfaces are non-compact. Therefore, in principle, the 
spatial integrals arising in this system diverge. In order to avoid this problem, one introduces a Euclidean 
fiducial 3-metric and restrict the integrals to a given cell, which is cubic with respect to the fiducial metric 
and whose fiducial volume is Vb fH*] . Owing to the imposition of homogeneity, the spatial diffeomorphism 
freedom is fixed and the associated constraints are trivially satisfied. Furthermore, the internal [SU(2)] 
gauge freedom is removed by choosing a diagonal gauge, in which the connection and the densitized triad 
are respectively given by 

K = cV-"'5l and = pV'^' . (1) 

So, one can identify the internal indices i = 1,2, 3 with those for the tangent space a — 1,2,3, and the 
three directions are equivalent owing to the isotropy. In this way, the two quantities c and p represent 
the only degrees of freedom of the geometry. This choice of parameterization for A]^ and leaves all 
the relevant physical quantities independent of the fiducial structures, as it is the case for the symplectic 
structure, determined by the Poisson brackets {c,p} = SttGj/S, where G is the Newton constant and 
7 > is the Immirzi parameter. 

In order to have non-trivial dynamics, it is necessary to add matter degrees of freedom. The simplest 
possibility is a homogeneous massless scalar field 0, with conjugate momentum p^, such that {(t>,p,p} = 1. 

The model is then subject only to the Hamiltonian constraint. Making use of the spatial homogeneity, 
the integrated form of this constraint for any lapse function N is C{N) = NC, where the Hamiltonian 
constraint 

C=-^-'VW\+SnG^ (2) 

has the standard densitization used in LQG. It is easy to check that, classically, both p and (f> are 
monotonous functions of the proper time. Usually, (f) is regarded as an internal time, with all quantities 
evolving as functions of it. Then, the classical solutions represent expanding or contracting universes, 
and all of them are singular. 



B. Quantum Representation 

In LQG, the gauge invariant information about the phase space is captured in holonomies of su(2)- 
connections and fluxes through surfaces. Similarly, in LQC one adopts as basic variables holonomies 

1/3 

along straight edges of oriented coordinate length ^iVq in the fiducial directions -that are given by 
hf{c) = e^'^'^^ for the direction i [l^]- and fluxes through squares normal to these fiducial directions 
-which are equal to p up to unimportant factors. The configuration algebra Cyl5 is the algebra of almost 
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periodic functions of c [2^, [23| generated by the matrix elements of the holonomies, A/"^ = gV=/2 ^^^^ 
momentum representation, the states defined by these matrix elements are denoted by The Cauchy 
completion of Cylg with respect to the discrete norm = 6^^^^' provides the gravitational part of the 

kinematical Hilbert space The operator p has a diagonal action on the basis formed by the states 

1^), whereas A/"^ shifts the state |//') to \n' + jj). 

The standard procedure in LQC to define the curvature of the su(2)-connection is to express it in terms 
of holonomies along closed loops and shrink their area to the minimum non-zero value A allowed by LQG, 
which is proportional to the Immirzi parameter and to the square of the Planck length Zpi = V Gh. Then, 

1 /3 

the condition that the physical area of a square with sides of (minimum) fiducial length /xV^ ' equals A 
Q leads to the operator relation 



1 



(3) 



Since the corresponding shift produced by Mfi is not constant in the basis of states |^), it is convenient 
to relabel these states by introducing an affine parameter which is proportional to the respective 

eigenvalue of the physical volume operator. In this manner, produces a constant increment in the 
volume, Sfp\v) = 



12/31 



+ 1), while 'p\v) = sign(u)(27r7;|,[vA|w|) 
On the other hand, for the matter part of the kinematical Hilbert space one takes the standard rep- 
resentation space of square integrable functions with respect to the Lebesgue measure, namely T-L^m ~ 



L (R, dcf). In total, the kinematical Hilbert space is "Hkin ~ H 



grav 
kin 



III. THE HAMILTONIAN CONSTRAINT OPERATOR 



Once the kinematical structure of the quantization has been introduced, we can proceed to represent 
the Hamiltonian constraint C as an operator. We will choose a symmetric ordering for the constraint 
which diff'ers from the orderings considered in Refs. 0, ll| but, as we will see, results more suitable. 

As usual in LQC, to represent the factor in the second term of Eq. ^ one appeals to Thiemann's 

procedure rewriting it in terms of Poisson brackets of holonomies with the physical volume :— \p\^^'^. 
On the other hand, in the gravitational part of the constraint, one expresses the curvature of the su(2)- 
connection using holonomies as explained above, in order to obtain a well defined operator 0-0] • 

In earlier works, the quantum Hamiltonian constraint was constructed along these lines, adopting cer- 
tain choices of factor ordering directly for the isotropic model In contrast, in our case, we will adopt 
an ordering motivated from previous studies performed in the anisotropic Bianchi I model [l7l [. In such 
a system, the sign of the triad plays an important role because there are three diflterent directions, rather 
than one. Then, products of two signs are not necessarily equal to the unity. As a consequence, in that 
case, one has to take carefully into account the presence of the signs when symmetrizing the constraint. 
To select a similar ordering in the present isotropic case, we start with the symmetric Hamiltonian con- 
straint constructed in Ref. jl7j and identify the three spatial directions. In this manner we arrive to the 
following constraint operator 



C 



1/2 
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^2 



+ SnGp; 



1/2 



(4) 



Here, the inverse volume operator is given by 



where 



\p\ 47r7?2yA 



T 






V 




_VW\_ 









(5) 



(6) 



The inverse volume acts diagonally on the basis states \v) and annihilates the state \v = 0). 
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On the other hand the operator 



n 



-I -1/2 



(7) 



coincides (in terms of holonomies and triad operators) with the one studied in Ref. , denoted there by 
Qi. The factors involving powers of the triad at the beginning and the end of this expression appear in 



a very particular manner, in analogy to the Bianchi I model [17|. Nonetheless, this concrete ordering is 
not relevant inasmuch as physical results do not depend appreciably on it, and one could select any other 
one. For instance, one could even ignore the quantum effects coming from the inverse of the volume [25j |. 
as was done in Ref. [sj . The really relevant features of our ordering, as we will see in the next section, 
come from the particular treatment of the sign of p that we have considered. 

It is straightforward to see that our Hamiltonian constraint operator annihilates the state |t; = 0) and 
leaves invariant its orthogonal complement, which we will denote by Ji^^^ . Hence, when studying the 
non-trivial solutions of the Hamiltonian constraint, we can restrict ourselves to this latter subspace. Note 
that "H^l^^ is just the Cauchy completion with the discrete norm of Cylg. the linear span of all the \v) 
states with v ^ 0. As in similar situations studied by us in previous works [T3, [H, HO, Hj] , the big bang is 
then resolved in the sense that the quantum equivalent to the classical singularity (namely, the eigenstate 
of vanishing physical volume) has been entirely removed from the kinematical Hilbert space. 

Once we have removed the kernel of the inverse volume operator, we can introduce the densitized 
version of C. Let us remember before that the physical states, which are annihilated by the Hamiltonian 
constraint, are generally not normalizable in the kinematical Hilbert space. We should seek them in a 
larger space. A natural home, as far as the gravitational part of the system is concerned, is the algebraic 
dual Cyl^ of the dense set Cyl5. In this dual space (tensor product any suitable space for the matter 
degrees of freedom), we can establish a one-to-one relation between any element ("01 annihilated by the 

(adjoint of the) operator C and any other element {ip'\ = {ipWl/V]^^"^ annihilated by the (adjoint of the) 
densitized version of the constraint, which is given by 



c ^-—n^ + 87tGpI. 



(8) 



This equivalent form of the Hamiltonian constraint is easier to impose since obviously and are 
Dirac observables which commute. 

Another way to get an easily solvable Hamiltonian constraint is to directly promote the classical 
densitized constraint C := VC to an operator. Nevertheless, as we have commented, such an object does 
not arise from the standard densitization of the constraint in LQG. For this reason, it seems natural to 
respect this latter densitization and show that one can construct a bijection between the solutions to the 
two considered Hamiltonian constraints. We emphasize that this bijection cannot be established in the 
kinematical Hilbert space, both because the physical volume operator is unbounded, and therefore it is 
not defined in the whole space, and because solutions do not belong to this space indeed. 



IV. CHARACTERIZATION OF THE GRAVITATIONAL CONSTRAINT OPERATOR 

As we have already seen, the operator f2^, which provides the gravitational part of the Hamiltonian 
constraint ([8]) , is a Dirac observable. Since the matter part of this Hamiltonian constraint is well known, 
the study of the properties of fl^, and in particular its spectral analysis, is the essential step for the 
resolution of the constraint. In this section, we will carry out a detailed analysis of this operator. We 
will identify first the superselection sectors which arise in the theory. Afterwards, we will perform the 
spectral analysis of fl^ in those sectors and determine its eigenfunctions. Remarkably, we can obtain their 
explicit expression, what will allow us to study their behavior. On the one hand, we will discuss how 
these eigenfunctions realize the commented no-boundary description. On the other hand, we will relate 
them with the eigenfunctions of the geometrodynamical (WDW) counterpart of fl^ for large v, in order 
to understand the WDW limit of our theory. Both features will give us insights about the existence of a 
quantum bounce. As a complement, we will also relate fi^ with the operator 17. Finally, we will compare 
our results with those of previous works on the polymeric quantization of the model, to point out the 
goodness of our prescription. 
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A. Superselection Sectors 



In LQC, owing to the discreteness of the volume representation, the gravitational part of the Hamil- 
tonian constraint turns out to be a difference operator instead of a differential one, as happens to be the 
case in the standard WDW theory. In the particular case of isotropic cosmologies, this operator produces 
a shift of four units in the label of the basis states \v). As a consequence, only basis states with support 
in discrete lattices of step four are related under its action. In addition, the Hilbert spaces of states which 
have support in such lattices turn out to be superselected, inasmuch as they are also preserved under the 
action of the complete Hamiltonian constraint and the physical observables. Whereas these results are 
general within isotropic LQC, in this subsection we will see that the lattices associated with our operator 
O^, and the corresponding superselection sectors, are simpler than those obtained in previous studies of 
the model. ^ 

The action of ft^ on the basis states 



of H^i^ takes the form 



h'\v} = -f+{v)f+{v + 2)\v + 4) + [fliv] + f_{v)] \v) - f-{v)f.{v ^2)\v-i), 



where 



f±iv) 



PI . 



iv±2)s±{v)g{v), 



s±{v) := sign(w ± 2) + sign(u). 



(9) 

(10) 
(11) 



and 




1 I 3 



- 1 



1 I 3 



if WT^O, 
if w = 0. 



(12) 



Notice that the combination of signs in the functions /±(w), which comes from the factor ordering that 
we have chosen to symmetrize the constraint operator, is the responsible of a remarkable property, namely, 
f-{v)f-{v - 2) = if u € (0,4], while f+{v)f+(v + 2) = for w e [-4,0). Therefore, the positive and 
negative semiaxes, u > and v are decoupled under the action of our operator fi^, as one can see from 
Eq. dni). In conclusion, this operator relates only basis states \v) with v belonging to one of the semilattices 
of step four £^ := {v — ±(£ + 4n), n £ N}, where e € (0,4]. In other words, fi^ is well defined in any 
of the Hilbert spaces Hf obtained as the closure of the respective domains Cyl^ := span{|w), v G Cf} 
with respect to the discrete inner product. Note that the non-separable kinematical Hilbert space HF^^ 



kin 



can be written as a direct sum of separable subspaces in the form V.^^^^ — (BiiHf < 



The Hilbert spaces Hf (8i L'^{M.,d4>) are preserved by the action of the Hamiltonian constraint (and 
then, as we will see, of physical observables). Thus, these Hilbert spaces provide superselection sectors 
in our theory. We can restrict the study to any of them. For concreteness, in the following we will work 
in ^^^^^(R, 

The difference with respect to previous works on the polymeric quantization of this model [lOl [Tl| 
is that the support of our sectors is contained in a single semiaxis of the real line, while in those works the 
sectors had contributions both from the positive and the negative semiaxes, which were not decoupled. 
Later on, we will discuss the advantages of working with our simpler sectors. 



B. Spectral Analysis 



Let us begin by showing that the symmetric operator (with domain Cy\f) is essentially self-adjoint. 
Actually, one can calculate that the difference between aJl^, with a :— 3/ {Air'y'^lp^h), and the gravitational 
constraint operator iJ^PS ol ^^f- CSl' both being defined in the Hilbert space iif (B'H^_g for e 7^ 4 (with 
natural domain Cyl-" U Cyl4_^), turns out to be a symmetric trace class operator. We obtain the same 
conclusion in the particular case e = 4 where, starting with the operator ail^ defined in © ^14, we 
have to define its action on \v = 0), e.g. by annihilation, since the operator -ff^ps does not decouple the 
state \v = 0). Taking into account that iJ^ps already proven to be essentially self-adjoint flo|, a well 
known theorem by Kato and Rellich ;26] ensures that so is fi^ as well. 
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In order to show that its restriction to, e.g., HJ is also essentially sell-adjoint, we will apply that, il 
A is a symmetric operator defined in certain Hilbert space H and p is any non-real number, then the 
operator is essentially sell-adjoint il and only if there exists no solution 10) € "H to the so-called deficiency 
index equation, A^cf)) — p\<j)) [27|. Let us suppose that il^ defined in Tif were not essentially sell-adjoint; 
this would mean that there exists a non-trivial solution to its deficiency index equation belonging to Hf, 
which in turn would provide a normalizable solution (identically vanishing in H^ -) when the operator is 
defined in the larger Hilbert space Tif ®'H^_g- We would then reach a contradiction because we already 

know that the operator is essentially sell-adjoint in this larger space. Therefore, 51^ has to be essentially 
sell- adjoint in T-L^ , as we wanted to prove. 

On the other hand, it was shown in Ref. [13] that the essential and the absolutely continuous spectra 
[27j ol the operator ^^aps both [0, oo). Once again, Kato's perturbation theory [2^ allows us to 
extend these results to our operator fi^ defined in 'H'^ ® ^4-£j since (up to a global lactor) it differs 
from H'j^pg in a symmetric trace class operator. In addition, taking into account the symmetry of il^ 
under a fiip of sign in v [f±{~v) ~ —f^{v)] and assuming the independence of the spectrum in the label 
e, we conclude that the operator D,^ defined in Hf is a positive (essentially) self-adjoint operator whose 
essential and absolutely continuous spectra are [0,oo) as well. Besides, as we will see in Subsec. IIVDI 
the (generalized) eigenfunctions of ft^ converge for large v to eigenfunctions of the WDW counterpart of 
the operator. This fact, together with the continuity of the spectrum in geometrodynamics, suffices to 
conclude that the discrete and singular spectra are empty psj . 



C. Generalized Eigenfunctions 



Let \e\) = J2vGcf ^A(^)k) denote a generalized eigenstate of fl"^ corresponding to the generalized 
eigenvalue A € [0,oo). For all n e N+, each coefficient e\{e + An) of this generalized eigenfunction turns 
out to be determined by the single initial datum e|(e) in the following manner 



ei(e + 4n) 



Gx{e-2) 



(13) 



where 



and 



-^W ■= -TTT' Gx{v}:=-- 



f+iv)' 



E 

0{a^b) 



Yl Fie + 2rp + 2) [] Gxie + 2s,) 



{rp} 



(14) 



(15) 



Here, 0{a — )■ b) denotes the set ol all possible ways to move from a to 6 by jumps ol one or two unit steps. 
For each element in 0{a 6), {rp} is the subset ol integers followed by a jump of two units, whereas 
{sq} is the subset of integers followed by a jump of only one unit. Note that F{i) — for all e < 2, so 
that in these cases the second term in Eq. does not contribute. 

As we stated above, the spectrum of is positive and absolutely continuous. In terms of a basis of 
generalized eigenstates \e\), the spectral resolution of the identity I in the kinematical Hilbert space 
is given by 



1 = 



dX\ei){e 



(16) 



Note that the integral runs just over the positive semiaxis and the spectrum is non-degenerate. The 
eigenlunctions satisfy the (5-normalization condition (e||e^,) = d{X — A'). This condition fixes the norm 
ol e| (e) in Eq. ((T3)) . The only remaining freedom is then the phase ol this initial datum. We finally fix 
this phase by taking e|(£) positive. The generalized eigenlunctions that form the basis are then real, a 

consequence ol the fact that the difference operator fl^ has real coefficients. 
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It is worth emphasizing that we have been able to solve the general eigenvalue equation of our grav- 
itational constraint operator 17^, determining explicitly the form of its generalized eigenfunctions. This 
contrasts with the level of resolution achieved in Ref. [7[ , where the generalized eigenfunctions of the cor- 
responding gravitational constraint operator were given in an iterative form and generated numerically. 
In this respect, we are in an optimal situation to progress in the comprehension of our system since we 
can now study its behavior analytically. In comparison with the exactly solvable model of Ref. Q , here 
we do not need to introduce simplifications in the system nor restrict the study to a particular sector 
of superselection, but our results are completely general. Remember that the construction of Ref. [1] 
was applied to a simplified version of the model and only in a specific sector whose support is centered 
symmetrically around v ~ 0. 



Another important issue that we want to investigate is the behavior of the quantum physical states 
(which we will determine in Sec. [V]) in the region of large volume. In particular, we want to discuss 
whether one recovers in that region the standard quantization performed in geometrodynamics, namely 
the WDW theory, whose predictions (for expectation values) in turn agree on semiclassical states with 
the classical ones obtained from General Relativity. To carry out such analysis, we will only need to know 
how the eigenfunctions of the gravitational constraint operator behave in the large v limit, since the 
other operator involved in the densitized Hamiltonian constraint, p^, has already been quantized in terms 
of the standard "Schroedinger-like" representation. In this subsection, we will obtain the eigenfunctions 
of the WDW analog of the operator il^ and relate them with the eigenfunctions e^^{v) for large v. 

As above for LQC, in the WDW quantization we work in the triad representation. The gravitational 
part of the kinematical Hilbert space of the WDW quantization can then be chosen as the space of square 
integrable functions of v with respect to the Lebesgue measure. The operator p acts by multiplication by 
the factorp = sign{v){2n^lp^^/A\v\)'^^^ , just as in the loop quantization, and the connection is represented 

by the derivative operator c — z 2(27r7Z|,j/A)-'^/^|u|-'^/^9t, so that [c,p] = ih{c,p}. 

Let us denote by Q_ the operator counterpart of the classical quantity (cp)^ in the WDW theory 
[defined in the Schwartz space 5(M)]. Since we want to compare its features with those of fi^, we choose 
for it the analog factor ordering, which gives rise to a symmetric operator that is well defined in the 
distributional sense: 



h = ~ 13^ ^\[sign{v)dy + 9„sign(w)] \v\ [sign(i;)a„ -I- a„sign(w)] = -^^ [l + 4vdy + 4(i;9„)2] , 



where /3 := Airjlp^. Note that we have simplified the expression of this operator by disregarding the 
non-contributing term |?;|(5(f ) in the second equality. 

^2 . 

Owing to well known properties of the operator —ivOy, we can ensure that fl is not only essentially 
self- adjoint in L^(R, dw), but also in each of the subspaces L'^{R^,dv). Hence, its action on the positive 
semiaxis v > and the negative one f < are decoupled, similarly to what happens with its analog 
in LQC, and we can restrict the study to L"^ {M.'^ , dv) . Furthermore, in this latter Hilbert space, the 

^2 

spectrum of [I is positive and absolutely continuous. Its generalized eigenfunctions, corresponding to 
any generalized eigenvalue A G [0, oo), can be labeled by a := ±\/A G M and are given by 



They provide a basis for L^(M+,(iw), normalized so that {e^\e^i) = i5(o' — a') (with S being the Dirac 
delta on the real line). 

^2 ^„ 
We see that the spectrum of fl has two-fold degeneracy, while we have shown that the spectrum of fl 

is non-degenerate. Therefore, any loop eigenfunction converges in the large v limit to a linear combination 

of the two corresponding WDW eigenfunctions (actually, one can rigorously prove that the limit of the 

loop eigenfunctions indeed exists [la])- Moreover, since the loop eigenfunctions are real, both WDW 

components must contribute with equal amplitude in that linear combination. Namely, the WDW limit 



D. Wheeler-DeWitt Limit 



(17) 




(18) 
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has the form 

el{v) -> r {exp [i(j)i{a)] e^{v) + exp [-i(t)g{a)] e_^{v)} , (19) 

where r is certain real number. In turn, one can check numerically that the phase shift 4>i{cr) has the 
following behavior 

0,(a)=r(|a|) + c, + i?,(|a|). (20) 

where T is a function of |cr| only, is a constant which depends on e, and linio— >.oo = (isj . 

So, whereas the eigenfunctions of our operator are determined by a single piece of initial data, they 
behave as eigenfunctions of a second order differential operator in the large v limit, therefore picking up 
a particular linear combination of the solutions to the eigenvalue problem of that differential operator. 
This nice feature of the polymeric quantization is the main responsible of the quantum bounce picture, 
together with the no-boundary description realized with our superselection sectors, as we will discuss in 
Sec. IVIl We postpone also to that section the comparison between the WDW limit of our theory and 
that of previous quantizations of the model. 

E. Operator il 

Once we have characterized the gravitational constraint operator f2^, let us relate it with fl to point 
out some interesting features which are due to the polymeric quantization and that do not appear in the 
analog WDW theory. We recall that (up to a multiplicative constant factor arising from a change in the 
basic commutators) Q coincides in fact with the operator 0^ extensively studied in Ref. :17]. We now 
summarize its properties. Like fl^, it is a difference operator, but with a constant step of two units in v 
instead of four. Its action is 

n\v) = -t[f+iv)\v + 2)- f^iv)\v ~ 2)] , (21) 

with f±{v) defined in Eq. ([10]). Taking into account that /-(w) = if w G (0,2] and f+{v) ~ when 
V £ [—2,0), we see that this operator does not mix states \v) with v belonging to different semilattices 
of the form '-^^Cf {v = ±{e + 2n), n G N}, where e G (0, 2]. In particular, our operator is well defined 
(with a natural choice of domain) in the Hilbert space '■'^^'Hf := T-C^^^ ©'H^^^^2- Furthermore, from the 
properties of O^, we infer that f] defined in ^^^Hf is an essentially self-adjoint operator whose spectrum 
is absolutely continuous, non-degenerate, and equal to the real line. In turn, its generalized eigenstates 
|e^), with support in '^'^^Cf and corresponding to the generalized eigenvalue a := ±-\/A G R, are formed 
by the direct sum of two generalized eigenstates of the squared operator fi^ for the eigenvalue A, one with 
support in the semilattice of step four C^^^ and the other supported in the semilattice ■ Explicitly 

one can see that \e%) = \/\^[\e^x} ® i sign(— (T)|e^"^^)] for cr 7^ 0, with (e^|e^,) = 5{(t — a') (like in the 
WDW case). For cr = 0, we define \e%^o) = \el^J) 

Let us comment that, whereas the eigenfunctions of fl^ have a well defined continuum limit for large v, 
those of r2 do not possess such a limit, since they are formed by two components, each of them admitting 
a WDW limit, but which are shifted by a phase equal to ±7r/2 (owing to the factor zLi in the linear 
combination). As a consequence, when v varies in '■-'^^Cf , the eigenfunctions e%{v) oscillate rapidly. This 
behavior is not present in the standard WDW theory, where the eigenfunctions of Q arc continuous in v 
for t; > [they coincide with those of £2 given in Eq. (IT51) ]. 

V. PHYSICAL HILBERT SPACE 

We can now complete the quantization by obtaining the physical Hilbert space and providing a complete 
set of observables. 

The matter term present in the densitized Hamiltonian constraint has been treated in a standard non- 
polymeric way. The (essentially) self-adjoint operator is positive with a two-fold degenerate spectrum. 
Its generalized eigenvalues are labeled by w^, with w G K. Let us call U the dense domain of definition of 
C, invariant under its action, from which one obtains the self-adjoint extension of this constraint operator 
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{U is the tensor product of, e.g., Cyl^ and a suitable domain for pp. Starting from this invariant domain, 

we can apply the group avera^ 
physical states have the form 



we can apply the group averaging method to find the physical Hilbert space 'H?^'^ [29l . I30j . The resulting 



^(«»=/ — TTei(«){V'+(A)exp[zc.(A),^]+V-(A)cxpH^(A)(/.]}, (22) 

where V'+ (A) and V'- (A) belong to the physical Hilbert space n^^y ^ L'^{M.+ ^uj-^{\)d\) and w(A) = 
[with a = 3/(47r7^Zpj?i)]. Regarding cj) as the internal time, we see that the solutions can be decomposed 
in positive (+) and negative (— ) frequency components, 5'±(w,(/'), which are determined by the initial 
data *±(v, 0o) through the unitary evolution 5'±(u, (j)) = U±{(j) — (f>o)'^±{v, 0o), where 

f/±((/)- 0o) = exp ±iVan^{(j)- 4>o)\. (23) 

A complete set of observables that allows us to interpret the system in an evolution picture is given by 
the constant of motion and the relational observable v\^^ (or if we do not restrict to w > 0, see 
Ref. Q). The latter measures the value of the volume when the time takes the value 0o- These Dirac 
observables preserve the positive and negative frequency sectors, so that, apart from the superselection 
already discussed, there exists further superselection with respect to the frequency. We can hence restrict 
the study, for instance, to the positive frequency sector. 



VI. DISCUSSION AND COMPARISON WITH PREVIOUS WORK 

In this section, we discuss the consequences of the quantization presented here as well as the main 
similarities and differences with respect to previous works on isotropic LQC, especially Ref. -7!], where the 
physical Hilbert space was originally determined and the quantum bounce of states that are semiclassical 
at late times was first studied satisfactorily. 

Let us remember that the analog of the operator il^ in Ref. 0], denoted by Q, is also (essentially) 
self-adjoint and positive. Furthermore, the quantum model constructed with O has superselection sectors 
supported in the lattices of step four C±\^\ := {v = ±|e| + 4n, n G Z}, where |e| € [0,2]. As we see, 
these lattices extend over the real line. In these superselection sectors, the spectrum of 6 is absolutely 
continuous and two-fold degenerate. Associated with this operator, a particular (5-)orthonormal basis of 
generalized eigenfunctions was chosen in Ref. Q- Their elements were denoted by efc(f) with fc G R, so 
that e\k\ (v) and e_|fc| (v) have the same eigenvalue. This basis is defined in such a way that e_|fe| {v) tends to 
S.-\k\i'") large positive v, where e_^,^^{v) is the generalized eigenfunction of the WDW operator analog 
of Q that provides an expanding wave. Similarly, we call e|j,|(w) the WDW eigenfunction corresponding 
to a contracting wave. The asymptotic behavior of e_|fc|(w) turns out to be given then by 

e-\k\ (v) ^ (v), e_|,| (v) (v) + Be^,^ (v). (24) 

Numerical analysis has shown that, for large \k\, A and B satisfy \A\'^ — \B\'^ = 1 and \A\ ~ \B\ ^ 1. 
As a result, the eigenfunctions e_\j^\{v) suffer an amplification in the negative semiaxis. In addition, this 
amplification is stronger as |fc| increases. 

On the other hand, since the parity transformation u — >■ — w is a (large) gauge symmetry of the theory, 
in Ref. [3] the analysis was reduced to the symmetric sector. Then, in general, it was necessary to 
join two different lattices so that the support of the states was symmetrically distributed around i; = 0. 
Thanks to the introduction of that symmetry, the analysis could be restricted to the positive u-semiaxis. 

Besides, in that work, the study was limited to the most interesting physical states: those which are 
semiclassical at late times. These states are provided by Gaussian profiles peaked at a large momentum 
of the scalar field or, equivalently, at a large negative k. Thus, in such regime the contribution of e\^{v) 
to the physical solutions is negligible. In this situation the form of the generalized eigenfunctions e\^{v) 
is irrelevant, and it was not necessary to calculate them. 

Once the parity symmetry is introduced and physical states are restricted in practice to the region of 
large negative fc's, one attains in Ref. a similar scenario to ours, in the sense that, for large positive 
V, the symmetric eigenfunctions that contribute significantly have a WDW limit in each lattice C±\^\ 
which is approximately of standing-wave type. This is just a consequence of Eq. (|24p and the commented 
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properties of the coefficients A and B, together with the implementation of the parity transformation. 
However, it is worth emphasizing that, while this standing-wave behavior is just an approximation vahd 
for fc <C — 1 in the case of Ref. , in our model this behavior is reached exactly and for all the eigenvalues 
of the scalar field momentum in the WDW limit. 

As for the procedure of Ref. to restrict to states in the parity symmetric sector, which leads to the 
mentioned union of two different lattices for generic |e| (namely, |e| ^ 0,2), this has some consequences 
which deserve special comment. Even when restricting the analysis to the sector k <^ —1, the WDW limit 
of the eigenfunctions carries a constant phase shift which depends on e. Actually, one can see numerically 
that the relative phase of the coefficients A and B presents the same kind of dependence found for our 
model in Eq. (I20[) (with e and |(t| now replaced with e and \k\). Therefore, even in the region of interest 
fc <C — 1, two different lattices possess different WDW limits and then their union does not admit a 
global limit. Remarkably, for the semiclassical states considered in Ref. 0, which are peaked for t; 1 
around two classical trajectories that do not overlap (one of an expanding universe and the other of a 
contracting one), the difference in the WDW limit of the two lattices is just a global phase for each of 
the two mentioned branches, which neither affects the norm of the state nor the expectation values of 
the observables. Nonetheless, this property of the WDW limit is not valid for more general states. 

The model studied in Ref. 0| was later simplified in Ref. mainly by disregarding the quantum 
corrections associated with the inverse volume operator. This simplification permitted to obtain an 
exactly solvable model in which the quantum bounce was shown in fact to be generic, although this 
result was attained, however, only for a specific superselection sector, namely, the one containing the 
state \v — 0). 

In comparison with Refs. 0,13)^ distinction of our proposal is that all the superselection sectors have 
a support contained in a single semiaxis. This allows us to restrict the study to e.g. v £ K+ in a natural 
way, without the need to appeal to a symmetrization process such as the one described above. In any case, 
note that the sectors of our model are perfectly compatible with the imposition of parity symmetry, which 
can be directly implemented by taking the direct sum of two sectors with support in the union of two 
semilattices, U Cj . Owing to the simplicity of our sectors, the spectrum of il^ is non-degenerate, what 
facilitates the exact and explicit calculation of the whole basis of generalized eigenfunctions [see Eq. (IT5|) ]. 
An analytical and numerical advantage of this non-degeneracy is that, to fix each eigenfunction, we only 
need to impose the positiveness of the initial datum e|(e), because its norm is completely determined by 
the normalization condition. As we have explained above, the WDW limit of these eigenfunctions takes 
the form ([T9l) . It is a combination of two WDW eigenfunctions which can be interpreted as contracting 
and expanding components (in v), or equivalently as incoming and outgoing components. They contribute 
with equal amplitude since the eigenfunctions are real, and in this sense the limit is exactly a standing 
wave. On the other hand, our eigenfunctions have support in a semiaxis which does not contain the 
potential singularity. This behavior does not arise from the imposition of any particular condition, like 
e.g. a boundary condition, but it is a natural feature of our model, explainable only by the functional 
properties of our gravitational constraint operator. From this perspective, we consider that our model 
provides an intrinsic no-boundary description. This implies that the outgoing component must evolve 
to an incoming component and vice versa, since the flux cannot scape across v = 0. Therefore, the 
expanding and contracting components must represent the two branches of a bouncing universe. Unlike 
in Ref. Q, where only a certain regime is considered, this result is independent of the spectral profile of 
interest. Furthermore, our analysis is valid for all choices of superselection sector (i.e., for all values of e), 
in contrast with the discussion carried out in Ref. i8]. In short, we obtain a completely generic quantum 
bounce. Obviously, the commented expanding and contracting components will be peaked (for large v) 
around well differentiated trajectories only for certain types of states. 

Let us conclude by remarking that the kind of no-boundary description that we have reached, which 
plays a key role in the arguments leading to the picture of a quantum bounce, is a characteristic of our 
model that is not shared by any of the previous works on fiat isotropic LQC [lOl. [Tl|. 

VII. CONCLUSION 

We have presented an alternative quantization in LQC of a flat, homogeneous, and isotropic universe 
in the presence of a massless scalar (homogeneous) fleld with the aim of improving our understanding 
of the theory. More explicitly, by a new choice of factor ordering, motivated by studies of anisotropic 
cosmologies, we have symmetrized the Hamiltonian constraint in a way which results specially appropriate 
to investigate some issues which had not been taken into account with sufficient care in previous analyses. 
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such us the precise relation between the Hamiltonian constraint and its densitized version, the possibihty 
to attain superselection sectors which, being as simple as possible, posses optimal properties, and a clear 
comprehension of the WDW limit of the theory. 

Our Hamiltonian constraint operator presents two nice features: under its action, the zero-volume 
state decouples, and states with different orientation of the triad are not mixed. Owing to the former of 
these properties, in our model the big bang singularity is kinematically resolved inasmuch as its quantum 
analog, namely the zero- volume state, is removed from the kinematical Hilbert space. This fact allows us 
to establish a bijection between the non-trivial solutions of the Hamiltonian constraint (with densitization 
equal to that used in LQG) and those of its densitized version (with respect to the volume). This bijection 
is not established in the kinematical Hilbert space, but in a larger space that provides a natural home for 
the solutions. For the densitized version of the constraint, the identification of some Dirac observables is 
straightforward. In summary, our analysis shows in a rigorous way how one can start with the density 
weight for the Hamiltonian constraint which arises naturally in LQG and make the passage to a simpler 
constraint with different densitization. 

The second feature commented above, namely the decoupling between triads with different orientation, 
has also important consequences. Thanks to it, our Hamiltonian constraint superselects sectors with 
support in a single semiaxis of the real line, instead of the whole real line. The simplicity of our sectors 
in turn simplifies considerably the construction of the physical Hilbert space in comparison with previous 
works [lol . [lH . In those references the superselection sectors have support in lattices of the real 
line. The study can be restricted to the positive semiaxis v > hy demanding parity symmetry, but in 
the generic case this is made at the cost of replacing the original individual lattices by the union of two 
lattices that are transformed one into each other under parity. In contrast, in our case the functional 
properties of the gravitational constraint operator restrict the study to t; > directly. As a consequence, 
the generalized eigenfunctions of this constraint operator, which are the elements that account for the 
quantum information about the geometry, provide a no-boundary description, arising in a single section 
of minimum non-zero volume v = s without the need to introduce any specific (boundary or large gauge) 
condition in order to affect the behavior in the vicinity of the origin. Another related consequence is that 
the spectrum of the gravitational constraint operator is non-degenerate, whereas the spectrum of the 
analog WDW operator is two-fold degenerate. Taking into account that the eigenfunctions of the loop 
operator can be chosen to be real, they must converge in the WDW limit (large v limit) to a standing wave 
composed of two WDW eigenfunctions equally contributing, one outgoing and the other incoming. In 
turn, the no-boundary behavior implies that neither the incoming flux nor the outgoing one can "scape" 
to the region v < 0. Thus, the only possibility is a bounce in which the incoming component becomes 
outgoing and vice versa. This shows the occurrence of a quantum bounce in a generic manner. Our 
conclusions do not only confirm the results obtained in previous works 0, Q j in particular concerning 
the robustness of the bounce, but also reinforce them inasmuch as the discussed scenario is completely 
general. We have neither focused the study on a concrete class of physical states, unlike in Ref. [3], nor 
simplified and particularized it to a specific superselection sector, unlike in Ref. Q. Finally, we note 
that the commented bounce scenario does not imply that the outgoing and ingoing components peak 
at non-overlapping classical trajectories, so that such a semiclassical behavior will be reached only for 
certain physical states (e.g. those considered in Ref. Q)- 

In conclusion, the simplicity of our quantum model has allowed us to solve it explicitly, gaining the 
physical intuition necessary to develop our understanding of LQG without resorting to numerical analyses. 
In this way, we have improved the control over the WDW limit and analyzed in depth the fundamental 
reasons behind the occurrence of a quantum bounce. 
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